Ab-initio determination of exchange integrals and Neel 
temperature in the chain cuprates 



A. B. van Oosten a and Frederic Mila 
a Max Planck Institut fur Physik komplexer Systeme, Nothnitzerstrasse 38, 01187 
Dresden, Germany 

Laboratoire de Physique Quantique, Universite Paul Sabatier, 31062 Toulouse France 

We report ab initio quantum chemical cluster calculations of the chain (J a ) 
and the largest interchain (J b ) Heisenberg exchange of the chain cuprates 
Ca 2 Cu0 3 and Sr 2 Cu0 3 . We find that J a is comparable to the in-plane J in 
layered cuprates and J a /J b ~250-400. Using recent theory we obtain close 
agreement with experiment for the staggered magnetic moments and 
critical temperatures. This implies that T N does not depend on the third 
parameter, J c «IJ b l, and cannot be calculated using spin-wave theory. We 
propose an explanation of this fact in terms of a 1D->2D cross-over. 



PACS numbers: 71.10.-w, 71.27.+a, 74.25.Jb, 75.10.Jm, 75.30.Et 



The one-dimensional (ID) spin 1/2 Heisenberg antiferromagnet is a nontrivial, but 
nonetheless integrable [1] quantum mechanical model with vanishing long-range 
magnetic order in the ground state. In real systems containing antiferromagnetic (AF) 
chains, interchain coupling often hinders the study of the ID behaviour. The linear 
cuprates Ca 2 CuC>3 and Sr 2 Cu0 3 [2] present a unique opportunity to compare exact 
results to experiment because of their exceptionally small ratio of inter- to intrachain 
coupling. The crystal structure consists of two interpenetrating, magnetically decoupled 
sublattices of infinite chains of S=l/2 [Cu0 3 ]" 4 units separated by Ca(Sr) +2 ions [3,4]. A 
single sublattice can be described by the Heisenberg hamiltonian 



where the summation runs over lattice vectors (ijk) and the unit cell convention is defined 
in Fig. 1. The magnetic susceptibility is entirely determined by the isolated chains, i.e. by 
the intrachain AF coupling J a . Two independent experiments yield incompatible values 
-J a «135-160 meV [3,5,6] and J a =-190 meV [7]. An even larger value of J a =-260 meV 
was estimated from the analysis of the mid-gap infrared absorption spectrum. For 
comparison, in the planar cuprates the in-plane coupling varies between -100 and -140 
meV[8,9]. A fascinating phenomenon is the transition to AF order at exceptionally low 
Neel temperature[10-12]. A recent theory developed by Schulz [13] predicts that the Neel 
temperature varies linearly with for isotropic coupling of the chains and that the square 
of the magnetic moment m % is proportional to the Neel temperature T N . This theory was 
applied to the problem of a two dimensional array of chains at zero temperature with the 
conclusion that such a system exhibits long range-order[14]. Recent muon spin rotation 
and neutron scattering experiments onCa 2 Cu03 and Sr 2 CuC>3 yield T N °c m n , 
approximately, and of the order of 1 meV. From a LDA band structure calculation J b 
has been estimated to be 3.6 meV and 0.8 meV for Ca 2 CuC>3 and Sr 2 Cu03, respectively 
[15]. However, Schulz's theory should not apply to the present case. It assumes 
isotropic interchain coupling, whereas in Ca 2 CuC>3 and Sr 2 CuC>3 the ferromagnetic [11] 
coupling along the c-direction should be extremely small, i.e. J c «IJ b l. Therefore, 
Schulz's theory predicts that T N and m K depend only on the average of J b and J c . This is 
in sharp contrast to the fact that, in the limit J c — > 0, m K -> at finite T according to the 

Hohenberg-Mermin- Wagner (HMW) theorem [16,17]. 

In this Letter we first determine the value of J a with an approach that is entirely 
independent from experiment and, for the quasi-2D cuprates, rivals with experiment in 
accuracy [8,9]. We then compute J b and show that Schulz's theory accurately predicts 
m n as well as T N , even though J c «J b . We then explain this unexpected result in a 

manner consistent with the HMW theorem. 

The computational approach for the intrachain coupling J a follows that of Refs. [8,9]. 
It is an accurately balanced calculation of singlet-triplet splitting of a cluster containing 



calculate the (MC)SCF singlet and triplet ground states of the planar, D 2 h symmetric, 
CU2O7 cluster that occurs in the chains. As demonstrated in Refs.[8,9], it is necessary 
and adequate to introduce local exchange and correlation effects on the bridging oxygen 
atom through subsequent admixture of excited (MC)SCF states, that differ from the 
MCSCF ground states by an O -> Cu electron excitation. We follow the non-orthogonal 
configuration interaction (NOCI) approach [19], which involves the computation of 
Hamilton and overlap matrix elements between determinants constructed from non- 
orthogonal orbital sets [20]. The NOCI approach leads to short, physically transparent 
wavefunctions and is size consistent. 

For the calculation of the interchain coupling J b , we perform all-electron (MC)SCF 
calculations of the singlet and triplet ground states of D 2 h symmetric C^Og and 
Cu20gM 16 (M=Ca,Sr) clusters (Fig. 2) consisting of two CUO4 units that are translated 
with respect to each other by one unit vector along the b-axis, so that they are located in 
next-nearest neighbour chains. Since there are no ligands, the exchange originates simply 
from the overlap between the Cu d-orbitals and can be found directly from the 
(MC)SCF ground state singlet-triplet splitting. For a study of the effect of the 
counterions, calculations were also performed on Cu20gM 16 (M=Ca,Sr). 

The clusters are embedded in a point charge environment with modified potentials at the 
nearest neighbour positions to the cluster. The complete specification of the background 
potential is available from the authors on request. For Cu and O the same contracted 
Gaussian basis sets are used as in ref. [9], namely a 14sllp6d -> 8s6p3d basis set for 
Cu and a 9s6p -> 3s3p one for O. It was found that J b is sensitive to extension of the Cu 
basis with diffuse d-functions, as will be discussed below in detail. For Ca we employed 
a 14s8p -> 5s2p basis set of double zeta quality [21] and a 12s6p -> 4s2p minimal basis 
set [22] and for Sr a (15s,9p,3d) -> (5s,3p,ld) minimal basis set [22]. 

Following Ref. [8,9], a first approximation to the intrachain coupling is obtained from 
a complete active space self consistent field (CASSCF) calculation of the singlet-triplet 
splitting of a D2h symmetric [CU2O7]" 10 cluster model. The active orbitals (whose 
occupation are allowed to vary) are even (g) and odd (u) linear combinations, d g and d u , 



of Wannier-like orbitals, d 1 and d 2 , localised at each Cu. The triplet wavefunction can be 
written as 



*F t = laadgdj = loodid 2 l . (2) 

Here a denotes the 0(2p a ) orbital at the bridging oxygen, which has the same symmetry 
as d u . The other closed shell orbitals are suppressed in the notation. The singlet 
corresponding to Eq.( 2) can be written as 



_ (1+S) lood g d g l-(l-S) laod u d u l 
V(2+2S 2 ) 



In Eq. (3) occurs the overlap S = <di I d 2 > as an additional variational parameter. As 
usual, the triplet and singlet SCF ground states are very well characterised by Cu +2 (3d 9 ) 
and O" 2 (2p 6 ) and the Cu holes have almost pure 3d(x 2 -y 2 ) character in a Mullrken orbital 
population. We find S-0.04 and E t -E s =20-30 meV. As before we admix to *P t and 
relaxed charge transfer (CT) excitations of the form 

¥ t * = ld u d u d g al , (4) 



, r J d g W- ld g d g d' u a'l 

V(2+2S* 2 ) 



In the excited singlet state ^ s * (5) d' u and a' are two non-orthogonal linear 
combinations of a and d u . As before d' u and o' are strongly overlapping, S*~0.55, and 
^s* describes a resonating covalent Cu-0 bond. The subsequent NOCI calculation 
involves a 2x2 non-orthogonal diagonalisation for the triplet and a 5x5 one for the 
singlet, since in the latter S and the coefficients defining d' u and a', namely S* and a 
d u /o mixing angle, are reoptimised. The resulting ground and excited states are separated 
by about 10 eV. For both Ca 2 CuC>3 and Sr 2 Cu03 a value of J a =— 1 19 meV is obtained, 



nearest Ca neighbours are included, as all electron atoms, the value for Ca2Cu03 

increases to 136.5 meV (Table I). From this it is plausible that also the Sr 2 Cu03 value 

will be between -119 and -136 meV. Our results corroborate analyses [5,6] of the 

susceptibility data [3] and are evidence against much larger values such as J a =-190 meV 

[7] and J a =-260meV [23]. 

The interchain exchange J b equals the singlet-triplet splitting of the Cu 2 Og based 

clusters shown in Fig. 2. There are no bridging O so that the (MC)SCF calculation 

following Eqs. (1) and (2) is sufficient. In this calculation the orbital d u transforms as 

b lu . Since the Cu-Cu distance of 3.278 A (Ca 2 Cu0 3 ) and 3.494 A (Sr 2 Cu0 3 ) is rather 

larger than a Cu-Cu inter-atomic distance, we study the effect on J b of basis set extension 

1 ? 

with diffuse Cu d- and O p-functions for [Cu 2 Og]~ clusters. The effect of diffuse d- 
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functions saturates with addition of oc=0.05, 0.016 and 0.005 bohr functions, while an 

_9 

extra oc=0.04 bohr O p-function has negligible effect (Table II). The diffuse functions 
contribute little to the total energy. The interchain coupling saturates towards values of 
J b =0.33 meV (Ca 2 Cu03) and J b =0.23 meV (Sr 2 Cu03). The overlaps of the d-orbitals 
are ~ 4 %o. The calculation was repeated for Cu 2 0g(Ca,Sr) 16 clusters shown in Fig. 2, 
using the extended Cu basis. The effect of the counterions is an increase of IJ b l from 0.33 
meV to 0.52 meV for Ca 2 Cu03 and from 0.23 meV to 0.28 meV for Sr 2 Cu03- The same 
results were obtained for the double zeta and the minimal Ca basis set. It is a reasonable 
assumption that further extension of the cluster will have negligible effect. 

From the fact that only the diffuse d-functions produce the antiferromagnetic interaction 
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it is clear that J b originates from overlap of the diffuse tails of the Cu 3d(x -y ) orbitals. 
In the language of the Hubbard model, one may describe this by a model system with a 
ground and a charge transfer excited state, connected by a matrix element 2t and 
separated in energy by U-V in the limit t -> [15]. Here U is the on-site repulsion in the 
d-orbital and V the Coulomb attraction between the electron and the hole it leaves behind. 
The ground state can be identified with the singlet (ldid 2 l-ldjd 2 l)/V2 and the excited state 
with (Ididjl + ld 2 d 2 l)/V2. Solving the 2x2 Hamiltonian and comparing to Eq. (3) one 
obtains the admixture of the charge transfer state as S=-2t/(U-V) and the energy lowering 



the relaxed non-orthogonal singlet and the overlap S, one can therefore determine the 
parameters t and U-V (Table IV). As one would expect, t decreases and U-V increases 
with increasing Cu-Cu distance. The increase of U-V of 0.3-0.4 eV can be ascribed to a 
decrease of V by 0.3 eV, as expected for two pointcharges. This suggests that V«4 eV 
and U>23 eV. This is larger than the atomic estimate E(Cu 3+ )+E(Cu + )-2E(Cu ++ )=47 
eV. Note, however, that U relates to the atomic 3d orbital, while only the diffuse tail of 
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the 3d(x -y ) orbital is involved in the hopping. The values for U and t are larger than 
values adopted in an LDA band structure calculation [15] values of 3-3.5 eV. 
Nonetheless, the large value reported here makes sense in view of the good experimental 
agreement obtained for J b . 

In Refs. [13,14] a theory of weakly coupled AF chains is developed. It is based on the 
response of a chain to a staggered magnetic field, % 1 S ^ > (T)=0.32/T ln 1/2 (5.8J/T) [24]. The 

interchain coupling is treated in the random phase approximation (RPA). In the present, 
anisotropic case the resulting expressions for m K and T N are 

m^= 1.017 VTjXTJT (7) 

and 

T N = 1.28 IJJ ln 1/2 | 4^ | . (8) 

Eq. (7) should be applicable to the present anisotropic case by taking 
IJ_L l=(l Jbl+U c l)/2=IJbl/2, as we can put J c =0 As can be seen from Table III, the resulting 
values of m K are in excellent agreement with experiment [12]. For the reasons stated 

above, Eq. (8) should not be valid in the present case. Nonetheless, the approximate 
experimental relationship T^m^ encouraged us to apply Eq. (8) anyway. The results 

are in surprisingly accurate agreement with experiment, as can be seen from Table III. 
The situation is graphically displayed in Fig. 3. 

We suggest an explanation of unexpected validity of Eq. (8). Roughly speaking [25] 
the 3D staggered susceptibility in the RPA can be written 



where Zj is the number of neighbours in the i-direction. This is equivalent to the set of 
expressions 
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As argued above, % st (T) cannot diverge at finite T, but is limited by fluctuations not 
present in the RPA treatment, so that Eq. (10) is not correct. For the case J a =J b it is 
known that % st (T) diverges as T" exp(AJ/T), where A is a numerical constant [26] and 
similar behaviour is plausible in the anisotropic case. Instead of Eq. (10) we propose the 
behaviour that is graphically displayed in Fig. 4. If J c is identically zero, the staggered 
susceptibility changes from % st to the unknown form % st at T , where T is still 
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determined by the condition l=z b IJ b l % st (T ). For non-vanishing J c , besides the cross- 
over at T*, a Neel transition occurs at T N <T*, where T N is determined by the condition 

2D 

l=z c IJ c l % st (T N ) whereby T N is guaranteed to vanish with J c as required by HMW. This 
is true as long as IJ C I is smaller than a critical value J*, determined by the condition that 
T N =T . From the fact that % S V(T)»% S V(T) unless T»T , it is clear that J c «J b . When 
IJ C I>J* one has T N >T* but, since % 2 S ^(T) rapidly reduces to X^CO for T>T*, T N will 
remain very close to T*. Thus in this case the Neel transition takes places very near to T* 
as given by Schulz's theory. In Fig. 5 an example of Schulz's prediction for T* and the 
linear spin wave (LSW) estimate for T N are plotted against J c . For IJ C I<J C the LSW 
estimate is qualitatively correct, while for IJ C I>J* it fails to reproduce the pinning of T N to 
T*. In this regime Schulz's prediction of T*=T N is quantitatively correct. Monte Carlo 
calculations on a plane system of 10 chains each containing 100 spins and with a 



Finally, it should be possible to experimentally observe the change in critical behaviour 
when the dimensional cross-over transforms into the Neel transition. Moreover, the 
behaviour discussed here should be universal and should also occur in systems other 
than magnetic ones that exhibit strong anisotropy and a transition from ID to 3D 
correlation. 

In summary, the intra- and interchain magnetic coupling, J a and J b , in Ca 2 Cu0 3 and 
Sr 2 Cu03 were computed. We obtain values of J a =- 119 meV for the intrachain coupling 
in both Ca 2 CuC>3 and Sr 2 Cu03, respectively, with Cu 2 C>7 clusters. This value is close to 
the value previously obtained for La 2 Cu04- For the Cu 2 07Mj5 (M=Ca,Sr) clusters J a is 
enhanced to -136 meV, which agrees with several analyses of the magnetic susceptibility 
of Sr 2 CuC>3. For the interchain coupling we obtain J b =-0.52 meV and J b =-0.28 meV, 
respectively. With these values Schulz-'s theory [13,14] predicts magnetic moments and 
Neel temperatures in close agreement with experimental data. Thus the transition 
temperature is determined by Heisenberg coupling in 2D only. We propose an, in 
principle, universal mechanism which pins the Neel temperature is to the dimensional 
(1D->2D) cross-over temperature T*, which is determined only by J a and J b , if the third 
coupling IJ C I exceeds a certain critical value J c «IJ b l. We suggest that an experimental 
study of the critical behaviour of Ca 2 Cu03 and Sr 2 Cu03 should be rewarding. 

We thank Stephan-Ludwig Drechsler, Walter Stephan and Karlo Penc for stimulating 
discussions. 



Tables 



a b c 

meV meV meV 

Sr 2 Cu0 3 (Cu 2 7 ) -32.8 -118.7 -100..-260 1 

Ca 2 Cu0 3 (Cu 2 7 ) -36.5 -118.8 — 

Ca 2 CuQ 3 (Cu 2 7 Ca 16 ) -48.1 -136.5 



Table I. Calculated intrachain coupling J a of the linear cuprates Sr 2 Cu0 3 and Ca 2 Cu0 3 : 
a) Triplet-singlet splitting of Eqs. (2,3); b) Non-orthogonal CI involving (2,3) and (4,5); 
d) Experiment, ^efs. [3,5-7,23]. 



added diffuse exponent 
2p Cu 3d 


Jb(meV) 
Sr 2 Cu0 3 Ca 2 Cu0 3 




-0.167 


-0.209 


0.05 


-0.176 


-0.252 


0.05, 0.016 


-0.226 


-0.321 


0.05, 0.016, 0.005 


-0.228 


-0.330 


0.04 


-0.166 


-0.213 



Table II. Calculated interchain coupling J b of Sr 2 Cu0 3 and Ca 2 Cu0 3 from Cu 2 Og 
clusters, for various basis sets. 





Jb(meV) 


magnetic moment (|i B ) 


Neel temperature (K) 






theory a experiment 


theory a experiment 


Sr 2 Cu0 3 


-0.276 


0.065 0.06(1) 


5.58 5.41(1) 


Ca 2 Cu0 3 


-0.522 


0.089 0.09(1) 


10.12 10 



Table III. Calculated interchain coupling J b along with ensuing theoretical values for 
staggered magnetic moment and Neel temperature in Sr 2 Cu03 and Ca 2 CuC>3. J], is 
calculated using Cu 2 0gMj5 (M=Sr,Ca) clusters. Theoretical values were obtained from 
the theory described in Refs. [13,14]. Experimental values were taken from Ref. [12]. 





2t 


U-V 




meV 


eV 


Sr 2 Cu0 3 (Cu 2 8 ) 


-75 


19.9 


Ca 2 Cu0 3 (Cu 2 8 ) 


-87 


19.6 


Sr 2 Cu0 3 (Cu 2 8 Sr 16 ) 


-83 


19.2 


Ca 2 Cu0 3 (Cu 2 7 Ca 16 ) 


-109 


18.8 



Table IV. Calculated hopping parameter t and on-site repulsion minus Coulomb attraction 
U-V of the linear cuprates Sr 2 Cu0 3 and Ca 2 Cu0 3 . 



Figure captions 



FIG. 1. The M 2 Cu0 3 (M=Sr,Ca) unit cell. Cu0 4 'plaquettes' are comnpleted with O 
outside the unit cell. 

FIG. 2. The Cu 2 8 M 16 (M=Sr,Ca) cluster. 

FIG. 3. Staggered magnetic moment and Neel temperature versus interchain coupling 
IJ b l, with J kept constant at 136.5 meV The solid curves represent the mean field 
expressions. The points correspond to Table III. 

FIG. 4. Staggered magnetic susceptibility in system of weakly coupled AF chains. 
Shown are % ®, % 2 S ^' RPA (Eq. 10) and an estimate of % 2 s ^that interpolates between 
% 2D,RPA for T>T ^ and 1/T exp(AJ/X ) with A-0.25 and J=VjJ^, for T<T*. 



FIG. 5. Neel temperature T N as a function of J c for J b =0.1 (all in units of J a ) as predicted 
by Schulz's theory and by Linear Spin Wave theory. The intersection of the two curves 
determines the LSW prediction of J*. 
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